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We analyze the spontaneous emission of a two-level atom interacting with a special class of struc- 
tured reservoirs of field modes with band gap edge coinciding with the atomic transition frequency. 
The exact time evolution of the population of the excited level is evaluated analytically through 
series of Fox-// functions. Over estimated long time scales, inverse power law relaxations emerge, 
with powers decreasing continuously to 2 according to the choice of the special reservoir. No trap- 
ping of the population of the excited level emerges. The same results are recovered in presence of 
N — 1 atoms, each one in the ground state, described by the Dicke model. The power of the inverse 
power law decay results to be independent of N. A critical number iVi*' is evaluated, such that, for 
N > iv£°, the inverse power law decay vanishes. 
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I. INTRODUCTION 

The influence of an external environment on the spon- 
taneous emission of an excited atom has been widely 
studied in literature. Great interest is devoted to the 
interaction of an atom with a reservoir of field modes, 
provided by cavities [l|, or dielectric media exhibit- 
ing photonic band gap (PBG) According to the 
structure of the reservoir, the atomic decay can be either 
enhanced or suppressed. Important effects emerge if the 
atomic transition frequency is near the PBG edge. 

In Ref. @ the decay of an atom placed in a three 
dimensional periodic dielectric, supplying a PBG 0], is 
analyzed. Under certain physical conditions involving 
the detuning, a portion of the initial population of the 
excited level is trapped and the dynamics is described 
by the photon-atom bond dressed or quasi-dressed states 
0, Over long time scales, the spontaneous emission 
exhibits an oscillating behavior depending on the detun- 
ing between the atomic transition and the band edge 
frequencies. If the transition frequency lies near (also 
outside) the band edge frequency, the trapping of the 
population of the excited level reveals [f| . 

The spontaneous emission of a two-level atom (TLA) 
has been also studied in presence of N — 1 unexcited 
atoms interacting according to the Dicke model [9| . The 
atomic decay of the population of the excited level ex- 
hibits a behavior similar to the case of one single atom, 
with an additional time scale factor equal to N 2 / 3 for 
an isotropic band gap and to N or N 2 for an anisotropic 
two dimensional or three dimensional band edges, respec- 
tively. For a detailed analysis we refer to [6] 

An attractive model for the inhibition of the atomic 
emission has been developed in Refs. [l(| and [ill, where 
the reservoir is dynamically modified through the inter- 
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action with an oscillating external bath. Physically, this 
condition is realized by making the main cavity leak to an 
external cavity with a moving mirror. The shape of the 
reservoir is unchanged but the couplings vary on time. 

Structured reservoir described by Lorentzian type 
spectral densities interacting with a multilevel atom has 
been analyzed in Ref. [T2|, [H[ and The cru- 

cial assumption that the frequency range is (— oo,+oo) 
makes the exact dynamics characterized by the poles of 
the spectral density in the lower half plane and the cor- 
responding pseudomodes [l5j]. A proper reservoir de- 
scribed by a spectral density which is the difference of 
two Lorentzians, exhibits a PBG. If the TLA is resonant 
with the gap, the population of the excited level is per- 
manently trapped. 

Recently [3, [13: an engineering reservoir approach 
has been adopted in order to delay the decoherence pro- 
cess affecting a qubit due to the interaction with its ex- 
ternal environment. Special reservoirs of field modes, in- 
teracting with a qubit in rotating wave approximation 
[HI , are designed in such a way that the qubit transition 
frequency coincides with a band gap edge. Over esti- 
mated long time scales, the decoherence process results 
in inverse power law relaxations, with powers decreas- 
ing continuously to unity, according to the choice of the 
special reservoir. 

The first theoretical approach to particle emission was 
made by Gamow at the beginning of the last Century, in- 
terpreting the decay as the appearance of purely outgoing 
waves in the solution of the corresponding Schrodinger 
equation at large distance and relating the absence of 
incoming waves to complex energy eigenvalues. Thus, 
the processes of tunneling and decay consist in res- 
onant states as solutions of the Schrodinger equation 
with purely outgoing boundary conditions. Survival and 
nonescape probabilities for particle decay exhibit inverse 
power law behavior with power 3 over long time scales. 
A detailed report on the argument and the above results 
is provided in Ref. [19| . 

In this scenario we study the spontaneous emission of 
an excited TLA interacting with the structured reservoirs 
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introduced in Ref. [17(. We aim to describe analytically 
the exact dynamics of the atomic population of the ex- 
cited level and compare with the decay provided by a 
periodic dielectric [g|. The atomic decay in presence of 
N — 1 atoms in the ground state, interacting according 
to the Dicke model, will also be analyzed. 



II. THE MODEL 

Wc consider a TLA interacting with a distribution of 
field modes in rotatin g w ave approximation and with 
dipole interactions fl2l TL3l \l§L 12(1 l2lj . By considering 
a system of units where h = 1, the Hamiltonian of the 
whole system is H = Ha + He + Hj, 



H A = wo|l) o(l|, H E = y^fc b[b k , 

k=l 

oo 

H I = i^g k (bl ®|0) a „(l|-6fc 



1 a a 



fc=l 



The state kets | 0) a and |l) a represent the ground and 
the excited state of the atom, while bt and b k are the cre- 
ation and annihilation operators, respectively, acting on 
the Hilbert space of the A:-th boson and fulfilling the com- 



mutation rule 



Sk,k' for every k, k' = 1,2,3,.. .. 
The constants g k , representing the coupling between the 
atomic transition and the fc-th field mode, depend on the 
geometry and are considered to be real without loss of 
generality. 

Starting from the initial state of the total system 



|*(0)) = |1)„®|0)j 



(1) 



where | 0).e is the vacuum state of the environment, the 
exact time evolution is described by the form 

oo 

|¥(t)) = c(t)|l) Q ® \Q) E + d ^)\ °) ® \ k )E, 

k=l 

\k) E = bl\0) E , k=l,2,.... 

The dynamics is easily studied in the interaction picture 
with respect to H$ + He, 

\V(t))i = e t(Hs+HE)t \^(t)) 

OO 

= C(t)\l) <g> \o) E + Mt)\o) ® \k) B , 



fc=l 



where I is the imaginary unity, C(t) 



A fc (i) = e %Ukt d k (t) for every k 
is described by the system 

oo 

C(t) = -J29kA k (t), 

k=l 

A*(t) = g k C(t)e«»*-»°», 



c(t) and 



1,2,.... The dynamics 



-i(u) k —u )t 



(2) 

1,2,3..., (3) 



solved in Ref. [6j. The amplitude (1| ® B (0||* (*))/, la- 
beled as C(t), is driven by the following convoluted struc- 
ture equation: 



d(t) = -(f*c) (t), 



C(0) 



i. 



(4) 



where / is the two-point correlation function of the reser- 
voir of field modes, 

oo 

} {t - t') = ^g 2 k e^ k - ua ^ t - t ') . 

k=l 

For a continuous distribution of modes described by 
t)(uj), the correlation function is expressed through the 
spectral density function J (uj), 



fir) 



J{uj)e- l ^-^ T dLU, 



where J (uj) = r\ (uj) g (uj) and g (uj) is the frequency de- 
pendent coupling constant. The above model was stud- 
ied by Garraway [H, EH and the corresponding exact 
dynamics was analytically described for Lorentzian type 
distributions of field modes. 



III. THE EXACT DECAY 

We analyze the exact time evolution of the population 
of the excited atomic level in case the structured reservoir 
is described by the following class of continuous spectral 
densities: 

2A (uj - upf 6 ((uj - up) /ujq) 

J a(u) = — 7 2 , (5) 



w ) 



A/a 3 - a > 0, 



1 > a > 0. 



The special spectral densities exhibit a PBG edge coin- 
ciding with the atomic transition frequency, an absolute 
maximum M a at the frequency (l a , 

M a = J a ((l a ) =Aa a ' 2 a a - 2 (2 - a) 1 "" 72 , 
(l a = uj + aa 1/2 (2~a) 1/2 , 

and are piece- wise similar to those usually adopted, i.e. 
sub-ohmic at low frequencies, uj/ujq > 1, and inverse 
power law at high frequencies, to/ ljq ^s> 1, similar to the 
Lorentzian one, though with different power, 



J a (uj) ~ 2 Aj 'a 2 (uj — uJo) a , 
J a (uj) ~ 2Auj a - 2 , lu^ 



UJ 

-oo. 



J i 



The convoluted structure equation (j4j has been stud- 
ied in detail in Ref. [171 ] for correlation functions /, gen- 
erated by the spectral densities ([5]). The particular case 
a = 1/2, characterized by a simple dynamics, has been 
separately analyzed in Ref. 16]. For the sake of clarity 
we report the main results and refer to the mentioned 
literature for a detailed proof. 
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A. The particular case a = 1/2 

The dynamics described by Eq. ((4]), takes a simple 
form in the special case a = 1/2. The solution is 



c(t) = ^ r Y, X iR(xi) e^r(l/2, X ^) 
V 71 " 1=1 

where R{z) is a rational function, 

( a i/2 + ,z) (t a i/2 + z ) 



i?(z) 



2^ ((1 + 1) a + 3 aV2 z + 2 (1 - i) z 2 ) ' 



(6) 



(7) 



while the complex numbers Xi)X2,X3 ancl X4i are the 
roots, distinct for every positive value of both A and a, 
of the polynomial 



Q(z) =Ti^/2/liA + iaz 2 + (l + i)a 1/2 z z + z 4 



(8) 



The analytical expressions of the roots are given by Eqs. 
([AT]) . (|A2]> . (|A3|) and |A4| of Appendix H Finally, the 
time evolution of the population of the excited state, 
P(t), given by the term | C(t)\ 2 , reads 



P{t) 



E» R(xi) e^r(l/2, X 2 t) 



1=1 



(9) 



The decay, fulfilling the initial condition P(0) = 1, is de- 
scribed by a simple quadratic form of Incomplete Gamma 
functions. 



B. The case 1 > a > 

In case the reservoir of field modes is described by the 
class of spectral density (JS]), the population of the excited 
level is the square modulus of the function C a (t), solu- 
tion of Eq. (|4]), where / is the corresponding correlation 
function of the reservoir. Following the arguments of Ref. 
171, the searched solution, 



[-l) n zl z r l - k t 3n ~ ak 



x HH 



n=0 k=0 
2 



k\{n-k)\ 



Zit' 



(-n,l) 
(0, 1) , (ak - 3n, 2) 



a 2 t 2 H^ 2 



(-n,l) 

(0, 1) , (afc - 3rc - 2, 2) 



(10) 



results in a series of Fox iJ-functions, defined through a 
Mellin-Barnes type integral in the complex domain, 



Tjm,n 
H p,<i 



(ai,ai) 



i (^p? ^p) 



1 

2vri 



, n™ x r {b 3 + fa) n™ =1 r (1 - a, - a,*) *- 

X / — r-=r; : — r US, 



IC 1L l=n+ 



i r(a I +a I «)lH =ro+1 r(l-ft i - / 9 i «) 



under the conditions that the poles of the Gamma func- 
tions in the dominator, do not coincide. The empty 
products are interpreted as unity. The natural num- 
bers m,n,p,q fulfill the constraints: < m < q, < 
n < p, and oti,/3j € (0, +00) for every i — 1, •■• ,p 
and j = 1, • • • , g. For the sake of shortness, we refer to 
22] for details on the contour path C, the existence and 
the properties of the Fox i?-function. The Generalized 
Mittag-LefHer [23|,[24|], the Generalized Hyper geometric , 
the Wright [25| and the Meijer G- functions (26| are par- 
ticular cases of the Fox iJ-function, thus, the amplitude 
C a (t) can be expressed as a series of each Special func- 
tion mentioned above. The parameters involved in Eq. 
(fTU| are defined as follows: 

Z\ — nAa 01 ^ 1 sec (na/2) — a 2 , z = iirAa a cos (na/2) , 
z a = -2nrAe-™ a/2 esc (na) . 

Simplified forms are obtained in particular cases. 
For example, the condition A = A^*\ where A^ — 
a 3 ~" cos {na/2) /tt, corresponding to z\ — 0, gives a 
power series solution, 

v n I 1 \ n r,\ v k r, n - k +3n— afc f 

c(*> a) = J2 E ■ • — - • - ' 1 



n=0 fc=0 



k\ (n- k)\T(3n- ak + 1) 



2 T (3n — afc + 1) 2 
a T (3n - afc + 3) * 



1 > a > 0. 



(11) 



In case the parameter a takes rational values, p/q, where 
p and q are distinct prime numbers such that < p < q, 
the solution of Eq. (j4|) can be expressed as a modulation 
of exponential relaxations [27j , 

/>oo />oo 

C p/q (t)= / (ft? / d^ p/q (ri,0e-^, (12) 
Jo Jo 

where 

n "if , / > 

/=i /c=i 

x sin (r) i 1/q sin (tt/q)) e n(&-"»(*7«)e 1/,r ) . 

The rational functions bi ^ (2) read 

d k-ij dz k-x ^ z g _ fl ) ( z9 + a ) ( z _ C,)" 1 ' IQ 



bi. k (z) 



(mi-k)\(k-l)\ 



for every I = 1, . . . , n, and k = 1, . . . , m;. The complex 
numbers & j • • • > Cn are the roots of the polynomial 



Q(z) = z 3 « + z lZ i + z a Z p + z 



(13) 



and mi is the multiplicity of Q, for every / = 1, . . . ,n, 
which means Q(z) = II™ =1 (2 — C/)™' an d m i = 3q. 

The particular conditions a = 3/4 and ^4 = 
a 9 / 4 cos (37r/8) /7r give a vanishing value of 21 and the 
roots Ci ! • • • > Ci can t> e evaluated analytically from the so- 
lutions of a quartic equation. For the sake of shortness, 
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the corresponding expressions are omitted. The remain- 
ing cases of rational values of a, force to a numerical 
evaluation of the roots of Q(z), once the numerical val- 
ues of both a and A are fixed. The necessary analysis of 
the function C a (t) is now completed. 

Starting from the initial condition ([1]), the interaction 
with a a reservoir of field modes described by the spectral 
density J a (w), gives the following time evolution of the 
population of the excited level: 



P a (t) 



x H 



1,2 



EE 

n=0 k=Q 
Zit 2 



( i \n „k ~n — k ±Zn—ak 

\ 1 ) z a z 1 

k\{n-k)\ 



(-n,l) 

(0, 1) , (ak - 3rc, 2) 



aH 2 



Zit' 



(-n,l) 

(0,l),(afc-3n-2,2) 



(14) 



The decay fulfills the initial condition P a (0) = 1 and gets 
a simplified form in the particular cases reported above. 



IV. INVERSE POWER LAWS 

The theoretical analysis of the exact decay, performed 
in the Section IIII1 leads to concrete results. The special 
reservoir described by the spectral density J\ii (w) gives 
an exact decay described by a quadratic form of Incom- 
plete Gamma functions and the corresponding asymp- 
totic expansions [2t| identify a time scale t and a decay 
factor C, 



max 



{|Xi| 2 , Z = 1,2,3,4}, 



(15) 



X, 



1=1 



such that, over long time scales, t/r >■ 1, the popula- 
tion of the excited level decays according to the following 
asymptotic form: 



P{t) 



t 



(16) 



The reservoirs described by the spectral densities 
J a (u>), for every a € (0,1), provoke an atomic decay 
resulting in inverse power laws over evaluated long time 
scales r a . For tj r a 3> 1 the population of the excited 
level, P a (t) is described by the asymptotic form 

P a (t) ~ C, a t~ 2(1+a \ i^+oo, l>a>0, (17) 
where 

4 a 2 a 4 ^-") esc 2 (ira) sec 4 (?ra/2) 
~ n 2 A 2 (r(l -a)) 2 ' 

A simple choice of the time scale is 



max < 1 , 



1/3 



1/a 



,3 



(18) 



Notice that the choice of the time scale for inverse power 
laws is not unique: in general, the time scales r and Tyn 
are different. 

The results obtained for the case a = 1/2 are confirmed 
by the numerical data plotted of Fig. [TJ Over long time 
scales, the population of the excited atomic energy level 
decays with an inverse power law, of power decreasing 
continuously to 2, according to the choice of the special 



reservoir. 



V. N ATOMS 

We extend the study of the spontaneous emission of an 
excited TLA in presence of AT — 1 ones in the ground state 
described by the Dicke model [l|, Q and interacting with 
the continuous distributions of modes (|5|) . Following Ref. 
Q , the Hamiltonian of the whole system reads 



00 00 

k=l 



k=l 
N 



Jl . m — 



n=l 



0(n) (n)( m \i 



I, m = 0, 1, 



where 1 1)(„) and | 0)(„) are the excited and ground state 
of the n-th atom, respectively. Notice that h = 1 and 
that the couplings between each atomic transition and 
the field modes are purely imaginary, i.e., the constants 
(/it are real for every k — 1,2,.... 

We adopt the following notation: the kets | J, M) 
are the normalized eigenstates of the operator J3 = 
{J2.2 ~ J\,i) /2, corresponding to the quantum number 
M, and the operator J 2 = Jf + (J2,i Ji,2 + Ji,2 ^2,1) /2, 
corresponding to the quantum number J. The N atoms 
are initially in the superradiant states \J,M = 1 — J ) 
with only one atom versing in the excited states and the 
reservoir, unentangled from the atoms, is in the vacuum 
states. The time evolution is described by the forms: 

*jv(t)>i = C N {t) \J,M = 1 - J) ® \0) E 

OO 

+ J2 A »,k(t)e-« u *-"°» \J,M = -J) ® \k) E , 

k=l 



C N (t) = -VN J2 9k A N ,k(t) e-'to-"-)*, 
fe=i 

A N>k (t) = VNgkC'N^e 1 ^"^, fc = l,2,3,. 



The above system is solved in Ref. [6j and leads to the 
following convoluted structure equation: 

C N (t) = -(C N *f N )(t), CW(0) = 1, (19) 

where /n is the two-points correlation function of the 
reservoir, 

00 

fw (t - f) = N^g 2 e -«("*-"o)(*-t'). 
fc=i 
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For a continuous distribution of field modes, the corre- 
lation function is expressed through the spectral density 

J (w), 

f N ( T )=N J (w) e-^ w - w °) T cL;. 
Jo 

The corresponding dynamics is obtained by replacing A 
with N A in the corresponding expressions obtained in 
Sections Hn] and HVI 

For the special case a = 1/2, the time evolution of the 
population of the excited level is described through the 
roots xna, XN,2, Xn,3 and xna of the polynomial 



!N 



(z) =Tr v / 2/7iNA + iaz 2 + (l + i)a 1/2 z 3 + z 4 , (20) 



obtained by Eqs. (|Alj) . (|A2|) . (jA3|) and (|A4|) of Appendix 
[A] by replacing the parameter A with AW. The popula- 
tion of the excited state results to be 



atom. For a = 1/2, a time scale tm and a decay factor 
Cat emerge, 



t n = max < |xjv,i| 



I = 1,2 



,3,4}, 



(24) 



such that, over long time scales, t 3> tat, the population 
of the excited level is described by the asymptotic form 



P N (t)~( N t- a , 
I . 4 

^ n = 4^ I X] R ( xn ' ( x7v - ; ) 



(25) 



As regards the general case, 1 > a > 0, a iime scoZe 
reveals, 



t n a = max < 1, 



^at, 



1/3 



3 £a 



l/a 



.3 



z Ar,0 



. (26) 



P N (t) 



1 I 

7T I 



x r 



(l/2, (XNj) 



4 

j'=i 

2 



(21) 



and fulfills the initial condition Pjv(0) = 1. 

The general case where 1 > a > is treated by re- 
placing the parameters z$, Z\ and z Q with zat, 0, ^at, 1 and 
zat.q, respectively, defined by changing A with AN, 

zn, 1 = 7rA7Va Q_1 sec (wa/2) — a 2 , zat :0 = Nzq, 
ZN,a = Nz a . 

The time evolution of the population of the excited level, 
PN,a(t) = I C N ^ a (t)\ 2 , reads 



P N ,a(t) 



EE 



ZN,lt 



n=0 k=0 
2 



(-N) n z k z£- k t 3n - ak 
k\{n-k)\ 



(-n,l) 
(0, 1) , (ak - 3n, 2) 



- a 2 t 2 



(-71,1) 

(0,1), (ak-Sn- 2,2) 



(22) 



and fulfills the initial condition Pat, a(0) = 1. 

Simplified forms are obtained in particular cases. For 



example, the condition A 



4W - „3-c* 
A^ 



defined as follows: 



(/■' "7rA cos(7ra/2) /n, corresponding to zm,i 
following pow 

00 n 

4*>) = EE 



0, gives the following power series solution, 

{-N) n n\z k a z^- k t 3n - ak 



n=0 fc=0 



M (n - k)\T (3n - ak + 1) 



x 1 



2 r (3n — afc + 1) 2 
° r(3n- afc + 3) * 



1 > a > 0.(23) 



such that, for tj tn >01 ^> 1, the population of the excited 
level exhibits inverse power law behavior described by 
the asymptotic form 



PN, a (t)~(N,at- 2(1+a \ t^+OO, 



C 



AT. 



4 a 2 a A{1 ~ a ^ esc 2 (no) sec 4 (na/2) 
tt 2 A 2 N 2 (r(l-a)) 2 ' 



(27) 
1 > a > 0. 



Notice that for N — > +00, the time scale for in- 
verse power law tends to the asymptotic value r 00j a = 

max |l, J 3/ (sin (na/2) cos 2 (ira/2)) 1 1//a /aj, while the 

decoherence factor vanishes, Cat, a — ^ 0, which mean that 
the inverse power law behavior vanishes as the number 
of atoms increase. We define a natural number Na^ such 



that C 



a 



2(l+ct) 



1, 



A W = 



2 a a 3 Q esc (7ra) sec (na/2) 



tiAY(\ - a) 



1 > a > 0, 
(28) 

where the square brackets indicate the integer part. The 
natural number Na estimates the critical order of mag- 
nitude for the number of atoms such that inverse power 
laws appear if the order of magnitude of N is less or equal 
to and vanishes if N > n£\ For every fixed value 
of the parameter a, the critical value can be enlarged by 
increasing the ratio a 3 ~ a /A. 

Again, the result obtained for the case a = 1/2, are 
confirmed by the numerical data plotted of Fig. [2j Over 
long time scales, the population of the excited atomic 
energy level decays with an inverse power law, of power 
decreasing continuously to 2, according to the choice of 
the special reservoir. 



VI. CONCLUSIONS 



As regards the dynamics over long time scales, it is 
described by inverse power laws as well as the case of one 



We study the spontaneous emission of a TLA inter- 
acting with a class of structured reservoirs exhibiting 
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|ln(P(t))| 




FIG. 1: The time evolution of [ln(P(£))|, the absolute value 
of the natural logarithm of population of the excited level 
obtained from Eq. ((9]), of a TLA interacting with the distri- 
bution of field modes J^/i (ui), given by Eq. ([5]) for a = 1/2, 
for < t/r < 10. The curve 71 (thick solid line) is obtained 
for A = (a/1000) 5/2 , 72 (thick dashed line) for A = 100 a 5/2 , 
73 (thin dashed line) for A = a 5//2 , 74 (thin solid line) for 
A = (a/2) 5//2 . The time scales for inverse power laws, given 
by Eq. (TPo) . read r ~ 1.40/a for 74, r ~ 0.84/a for 73, 
r ~ 0.06/a for 72, r ~ 7121. 40/a for 71, respectively. 

|ln(P„(t))| 




a t 



20 40 60 80 100 



FIG. 2: The time evolution absolute value of the natural log- 
arithm of Pff(t), given by Eq. (|21[) . the population of the 
excited energy level of a TLA interacting with a continuous 
distribution of field modes described by J\i% (w), given by Eq. 
(5) for a = 1/2, in presence of N — 1 atoms in the ground 
level, described by the Dicke model, for < at < 100 and 
A = 10- 4 a 5/2 . The curve Ci (thin solid line) corresponds to 
N = 3, C 2 (thin dashed line) to N = 12, C 3 (thick dashed 
line) to N = 60 and C 4 (thick solid line) to N = 500. The 
corresponding time scales for inverse power laws, given by Eq. 
(26) for a = 1/2, read r 100 oo,i/2 - 0.84/a, r 60 ,i/ 2 ^ 46.1/a, 
Ti2,l/2 — 206. 9/a, — 789. 0/a. The critical number, 

given by Eq. (25) for a = 1/2, is AT}*' = 3591. 

a PBG edge coinciding with the atomic transition fre- 
quency, and piece-wise similar to those usually adopted, 



i.e., sub-ohmic at low frequencies and inverse power laws 
at high frequencies, similar to the Lorentzian one. The 
exact time evolution of the population of the excited level 
is evaluated analytically through series of Fox-H func- 
tions. Over evaluated long time scales the population 
of the excited level vanishes according to inverse power 
laws, with powers decreasing continuously to 2 depending 
on the choice of the special reservoir. 

We also consider also the decay of an excited atom in 
presence of N— 1 atoms, each one in the ground state, de- 
scribed by the Dicke model, interacting with the special 
class of reservoirs. Again, the exact time evolution of the 
population of the excited level is described analytically 
through series of Fox-H functions and inverse power laws 
emerge over estimated long time scales. Though the pow- 
ers are not affected by TV, the inverse power law decay 
tends to vanish if N ^> Na , where Na is a properly 
defined critical number. No trapping of the population 
reveals. 

The above results are of interest if compared to the fol- 
lowing ones. An atom interacting with a radiation field 
created by a three-dimensional periodic dielectric with 
transition frequency lying near the edge of the PBG, ex- 
hibits an oscillatory decay of the population of the ex- 
cited level with trapping, spectral splitting and subnat- 
ural linewidth. In case the atom is surrounded by N — 1 
atoms in ground states, the time scale factor of the decay 
is iV 2 / 3 for an isotropic band gap, N or N 2 for anisotropic 
two-dimensional or three dimensional edges, respectively 

The special class of reservoirs of modes can in principle 
be realized with materials providing the PBG structure. 
A A^-period one dimensional lattice provides a band gap 
through an appropriate sequence of dielectric unit cells 
[28| . The corresponding density of frequency modes de- 
pends on the transmission coefficients of each unit cell. 
Also, the advanced technologies concerning diffractive 
grating and photonic crystals allows the realization of 
tunable ID PBG microcavities [H-[3l|. The action of 
such a structured environment on a qubit could be a way 
of delaying the decoherence process with fundamental ap- 
plications to Quantum Information processing Technolo- 
gies. 



Appendix A: details 

The exact expressions of the roots of the fourth order 
polynomial Q(z), given by Eq. (fI15|) , are described by the 
following forms: 
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where the parameters 7, and A are defined as follows: 
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